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Abstract 

The quasi-species equation describes the evolution of the probability that a random individual 
in a population carries a given genome. Here we map the quasi-species equation for individuals 
of a self-reproducing population to an ensemble of scalar field elementary units undergoing a 
creation and annihilation process. In this mapping, the individuals of the population are mapped 
to field units and their genome to the field value. The selective pressure is mapped to an inverse 
temperature /3 of the system regulating the evolutionary dynamics of the fields. We show that the 
quasi-species equation if applied to an ensemble of field units gives in the small (3 limit can be put in 
relation with existing stochastic quantization approaches. The ensemble of field units described by 
the quasi-species equation relaxes to the fundamental state, describing an intrinsically dissipative 
dynamics. For a quadratic dispersion relation the mean energy (U) of the system changes as 
a function of the inverse temperature /?. For small values of /3 the average energy (U) takes a 
relativistic form, for large values of /3, the average energy {U} takes a classical form. 

PACS numbers: 



1 



I. INTRODUCTION 



Increasing evidence shows that the Darwinian theory of biological evolution de- 
scribed by the quasi-species equation shares surprising similarities with quantum me- 
chanics at the mathematical level 4hJl6l] . The quasi-species equation describes the evolution 
of self-replicating macromolecules such as RNA or DNA or asexual organisms. In order 
to solve the quasi-species equation different methods from quantum mechanics have been 



used, including quantum spin chains {4], path integrals |5|-|7| and the Schrodinger equation 
in imaginary time js-ll]. This intriguing relation extends also to the steady state of biologi- 
cal evolution which is described by a Bose-Einstein distribution in a number of evolutionary 



0£i 



models of sexual and asexual populations 



I2J-I14J . Creation-annihilation operators and Fock- 



space formulations have been used by several authors for modelling the stochastic dynamics 
of biological evolution 



15 



16] . Finally many-body theory approaches have been applied to 



the description of stochastic gene expression 



The relation between quantum mechanics and stochastic dynamics has been deeply explored 



over the years 



18] . In particular the existing theory of stochastic processes for diffusing par- 



ticles is described by Langevin and Fokker-Planck equations. The latter are known to be 
related to the solution of the Schrodinger equation in imaginary time. Thus the stochas- 
tic diffusion dynamics is used in the existing stochastic quantization approaches to define 
stochastic processes whose probability density converges to path integrals of some specific 



quantum system |19l-l21|. 

Biological evolution is essentially a stochastic process determined by the reproductive rate 
(depending on the Fisher fitness of the individual and on the selective pressure) and by 
the mutation rate. The resulting stochastic birth-death processes are described by the 
quasi-species equation which is the mean-field equation for this stochastic dynamics. The 
quasi-species equation for biological evolution is the simplest example for showing that many 
biological systems have an error threshold. If, in presence of a single-peak landscape, the 
mutation rate is higher than a threshold value the biological system is not able to climb 
to the fitness peak of the energy landscape. In addition to the mean-field description of 
biological evolution give n by the quasi-species equation many interesting effects are due to 



the stochastic dynamics 



22|, fluctuating environment j^J, unbounded fitness landscapes 23] 



or spatial structure and time delay effects 24]. 
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Already in 1975 J. B. Anderson 25] proposed to use a birth-death process, i.e. a Monte 
Carlo with cloning, to reproduce the ground state of a given quantum system and the meth- 



ods is currently used in the literature of optimization [26|, |27j. The starting point for our 
investigation is that the quasi-species equation defines stochastic processes that in the small 
selection limit can be related to stochastic quantization a la Anderson 25[. These classes 
of stochastic processes involve creation-annihilation processes and a dissipative dynamics 
captured by the quasi-species equation with different functions playing the role of the Fisher 
fitness and the mutation rate. 
In a previous paper 



111 ] we have investigated the mathematical structure of the quasi-species 
equation by applying it to an ensemble of particles undergoing a creation-annihilation pro- 
cess under the drive of an harmonic potential playing the role of the Fisher fitness and a 
Gaussian noise playing the role of the mutation rate. Moreover the inverse temperature 
(3 plays the role of the selective pressure. This process can be viewed as an equivalent of 
the Ornstein-Uhlenbeck process in which the particles instead of diffusing are undergoing 
a creation-annihilation process. In 11] we found that the probability distribution of the 
process can be decomposed into eigenfunctions corresponding to a a discrete spectrum of 
eigenvalues. This spectrum is equivalent to the spectrum of the quantum mechanical har- 
monic oscillator in the limit of small inverse temperature /3. 

In this paper we extend these results by investigating the evolution of an ensemble of ele- 
mentary units, that we call field units, associated with a scalar field p(x). As individuals 
in a population carry a genome {cr^ } i=li where Oi indicates the nucleotide at position i, 
the field units that we introduce in this paper, are associated with a field p(x). These field 
units might represent for example a strategy of an agent in an evolutionary game theory. 
The field units undergo a creation-annihilation process inspired by biological evolution. A 
precise mapping can be done between the evolution of field units and biological evolution 
of asexual populations. In this mapping each field unit corresponds to an individual of the 
population and each value of the field corresponds to the possible genomes of the individual 
in the biological population. Moreover the energy of the field unit corresponds to the Fisher 
fitness of biological evolution and a Gaussian noise corresponds to random mutations. Fi- 
nally the inverse temperature /3 corresponds to the selective pressure of biological evolution. 
Solving the quasi-species equation for such a model implies that we need to find an expres- 
sion for the probability P({p(x)}) that a field unit carries a field p(x). This quantity is a 
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functional of p(x) and therefore we need to extend the quasi-species in order to describe the 
time evolution of a functional. For simplicity we solve explicitly the model for the energy 
of a single field unit which is quadratic in the field strength, but the same framework might 
be extended to other form of the energy U. We find that the distribution of the field units 
can be decomposed in a discrete set of eigenvalues. These eigenvalues are characterized by a 
set of discrete numbers {A^} extending in this way the stochastic quantization of particles 



found in 



ll| to a stochastic quantization of fields. Interestingly we find that in the limit of 



small inverse temperature f3 the average energy (U) takes a relativistic form while for high 
value of the inverse temperature the average energy (U) takes a classical form. Finally the 
dynamics is an off-equilibrium dissipative dynamics that asymptotically in time converges to 
the non-degenerated fundamental state {N^ = 0}. This aspect of the evolutionary dynamics 
is especially interesting since it might allow application of this dynamics for devising new 
optimization algorithms. 



II. THE QUASI-SPECIES EQUATION OF BIOLOGICAL EVOLUTION 

The genome of an asexual organism is formed by a single copy of each chromosome. For 
a given population the genome is determined by the sequence {a} = (a±, cr 2 , . . . 0%, . . . a N ) 
where % = 1,2, ... ,N is the number of genetic loci, and each variable <7j can take four 
possible values <7, = 1,2,3,4 indicating the nucleotide at each genetic locus. According 
to Darwinian evolution a biological population evolves under the drive of selection and 
mutations. We associate with a genome {a} the Fisher fitness U({a}) and a reproductive 
rate W({cr}) = exp[— (3U({a})] where j3 indicates the selective pressure. The selective 
pressure f3 is such that for — all genotypes have the same reproductive rate W({<r}) = 1, 
and when (3 ^> 1 the different genotypes corresponding to different Fisher fitness U({a}) 
have very different reproductive rate. In order to simplify the model, we assume that the 
death rate is a random process independent of the genome of the individual. The quasi- 
species equation Q is the mean-field equation for the probability distribution P({a}, r) that 
at time r an individual has genome {a} and is given by 

dP({a},r) 



dr 



M 



W}\W'} 



P({a'},r) -Z T P({a},r) (1) 
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where Z T is a normalization constant and can be expressed as Z T 



-0U(W}) 



, i.e. as 



the average over the probability distribution P({a},r). The operator Mjo-mo-/} in Eq. (CQ) 
describes the occurrence of mutations and acts on a generic function h({a'}) as 



M W}lW} h({a'}) = £Q(MIK})MK}) 

{*>} 



(2) 



The matrix Q({o"}|{o''}) represents the probability of mutations. If /i is the mutation rate 
we have 



Q(WIK}) = I1 

i 

with <5(<7i,<7j') indicating the Kronecker delta. 



1 - iM)6(<T i7 afi + - 



(3) 



III. EVOLUTION OF AN ENSEMBLE OF PARTICLES 



In JjJ we applied the quasi-species equation to an ensemble of particles that are located 
in one-dimensional space and that undergo a creation-annihilation process. The evolution 
of the particles can be mapped to biological evolution. In this mapping, the fitness function 
corresponds to the energy of a particle, and mutations correspond to a stochastic noise. 
Finally we assume that the probability to find a particle at a certain position obeys the 
quasi-species equation. 

If P(z, t) is the probability that a particle is at position z at time r, the quasi-species 
equation for this ensemble of particles reads 

dP(z,r 



di 



M, 



e -W)p(V )r ) -Z T P(Z,T 



(4) 



where the partition function Z T is given by 



dz' [ dzQ(z,z')e-P u( - z,)) P(z') 



(5) 



and the operator M Zi2 ,', applied to a function e f(z), acts as 



M 2 , 2l 



dz'Q(z,z')e-^ z,) f(z') 



(6) 



where Q(z, z') describes the stochastic noise playing the role of mutations for the evolution 
of the ensemble of particles. The inverse temperature (3, with /3 > 0, plays in this stochastic 
dynamics the same role as the selective pressure in the biological evolutionary dynamics. 
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For simplicity we assume that the position z at time t + dt is related to the position z' at 
time t by z = z' + rj where t] is a noise with Gaussian distribution. Therefore we take 



QO, z') = ^ ^ y dr)5(z' -z- r))e 



(7) 



where 0/D is the variance of the noise. In ll| we studied this process and we found that the 
probability of a certain configuration of the ensemble of particles can be decomposed into 
eigenf unctions of the evolutionary operator corresponding to a discrete spectrum of eigen- 
values. For low values of /3, the quasi-species equation is directly related to the Schrodinger 
equation in imaginary time. In fact, introducing in Eq. (112jl the quadratic expression (ED) 
and using the Fourier representation for the delta function we get 



M, 



(8) 



where M = ^J yTsd an( ^ ^he Hamiltonian H(x, k) of this system is identified as H(x, k) 



y/c 2 + U(x). The action of the operator Mj^ on a function f(x) for (3 1 is given by 



-0U(x') 



M, 



Therefore for (3 1 the evolution Eq. ( TTOT) can be written as 



~J\f dx' dk{l- (3H(x, k)}e 



ik(x-x') ft I 



fix 1 ) . 



dP(x,t) 
~dt 



P 



Dd 2 P(x,t) 

y dx 2 



U(x)P(x,t) 



+ (^-l)P(x,t) . 



(9) 



This equation is particularly interesting since it is linear in the Hamiltonian H (x, k) and 
can be put in relation with the Schrodinger equation in imaginary time. Therefore in this 
limit we are considering a quantum mechanical system evolving in imaginary time and we 
recover the classic model introduced by J. B. Anderson 25| in order to find the ground state 
of a quantum mechanical system. 

Nevertheless in the high-energy limit /? ^> 1 the quasi-species equation strongly deviates 
from the Schrodinger equation. The case of an harmonic potential has been explicitly solved 



in 



11] where it was shown that the eigenvalues of the evolutionary operator in the |3 < 1 



limit are given by the spectrum of the quantum harmonic oscillator. 
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IV. EVOLUTION OF FIELD UNITS 



In this paper we apply the quasi-species equation to a discrete ensemble of field units 
of a scalar field. These field units generalize the concept of an individual of a population 
to a wider context. As the genome of an individual is defined by a sequence of variables 
{a} we assume that each field unit is associated to a value of the scalar field p(x) taking 
complex values and defined on a continuous one-dimensional variable x. As the individuals 
in a biological population undergo a birth-death process we assume that the field units 
undergo a creation-annihilation process and we suppose that they follow the quasi-species 
equation. In the mapping between the population of asexual individuals and the ensemble 
of field units, the Fisher fitness is mapped to the energy U({p(x)}) associated to the value 
of the field, and the mutation probability is substituted with a Gaussian noise. We might 
describe this dynamics in the framework of an evolutionary game context. We assume that 
the elementary units are individuals with a given strategy p(x) and they play against a 
common player with fixed strategy playing the role of a fixed environment. The probability 
that a strategy p(x) reproduces depends on its fitness (payoff of the strategy in the given 
environment) W ({p(x)}) = e~P u " p ( x "> and there is a random drift of the elementary units 
indicating that strategy might fade away if not successful. We indicate by the functional 
P({p(x)},r) the probability that at time r a field unit is associated with the field p(x). 
The quasi-species equation for the ensemble of field units is an equation for the functional 

p ({p( x )}) S iven by 



dP({p(x)},r) 



M {p(x)}i{ p (x)} [e 
Z T P({p(x)},r) 



•^««*)})p({ / 5(x)},r) 



dr 



(10) 



where the partition function Z T is 



J Vp(x)Vp*(x)Vp(x)Vp*(x) x 
xQ({p(x)}, {p(x)})e-? u «W»P({p(x)}, r). 



(11) 



The operator ^■{ P (x)},{p(x)} applied to a function h(p(x)) acts according to 




xQ({p(x)},{p(x)})h({p(x)}) 



(12) 



7 



where the kernel Q({p(x)}, {p(x)}) describes the stochastic noise playing the role of muta- 
tions. The inverse temperature (3 > plays in this stochastic dynamics the same role as the 
selective pressure in the biological evolutionary dynamics. For simplicity we assume that the 
function Q({p(x)}, {p{x)}) has the form of a noise with Gaussian distribution. We assume 
that the field p(x) is defined on a box of size L with periodic boundary conditions. Therefore 
we take 

Q({p(x)}, {p(x)}) = Afe~^ I dx[p*(x)- P *(x)][p(x)- P (x)] ( 13 ) 

where M is the normalization factor. The variance of the noise Q({p(x)}, {p(x)}), given by 
1/C = j3mc 2 , is proportional to the inverse temperature and to a typical energetic scale of 



the system that we indicate by mc 2 28(. The energy U({p(x)}) defines the probability that 
one field unit replicates and generates another field unit. The lower the energy U({p(x)}) 
the higher is the probability that a field unit replicates. By repeating the same argument 
used in the preceding paragraph, the equation (TTOT) reduces to a Schrodiger equation for 
the field functionals only in the limit of small values of 0. For large value of the selection 
pressure the dynamical equation we are considering strongly deviates from it. 

In order to derive a specific solution of the stochastic dynamics of the ensemble of scalar 
field units described by Eq. ([TO]) , we take the energy quadratic in the field p(x). In particular 
we take 



h 2 l 

U{{p{x)\) = 

with 



U({p(x)}) = --i- / dxC(p(x)), (14) 



The choice of this energy might be generalized to energy including higher powers of the field 
p(x) or to consider adaptive fitness landscapes inspired by models of pattern formation or 



models of motion of active walker driven by chemotaxis [29] . Since p(x) is defined on a box 
of size L, it can be decomposed into a discrete series of Fourier components p(k) associated 
with a k = that can be decomposed according to p{k) = pB,{k)+ipi{k) where pii(k), pi(k) 
are the real and the imaginary parts of p(k). The energy U({p(x)}) = U({p(k)}) given by 
Eqs. (Tl4l) . ([15]) can be written as 



U({p(k)}) = -L. £ e 2 (k)\p\{k) + p){k)] (16) 

where e{k) is defined as 

e 2 {k) = [(hkc) 2 + (mc 2 ) 2 }. (17) 



We note that we can express also the noise Q({p(x)}, {p(x)}) = Q({p(k)}, {p(k)}), given by 
Eq. ( 1T3|) . in terms of pn(k), pi(k), obtaining 

Q({p(k)},{p(k)})=U II J^e- C ^-^\ (18) 

k a=R,I ' 

In order to solve the evolution of this ensemble of field units given by Eq. (flOl) . we study 
the eigenfunction problem 



M 



{p(x)},{p(x)} 



e 



m{ ^ x)}) nn(p(x))} = \rt*n({p(?)}) (19) 



where n = {nji(ki), rijiki), . . . , n^(/ci), n/(/c £ )} characterize the eigenfunctions which are so- 
lutions of the eigenvalue problem in Eq. f l2Tj) . The eigenfunction 7r^({p(A;)}) of the eigenvalue 
problem Eq. (ED), with the choice of U({p{k)}) given by Eq. ( fl6l) and Q({p(k)}, {p{k)}) 
given by Eq. ( TTBl) factorizes on functions defined on the single variables pnj(k), i.e. 

Kn({p(k)}) = U II ^n a (k)(p a (k)) . (20) 

k a=R,I 

The single functions ^n a (k)(Pa(k)) with a = R,I are given by 

7r na{k) (p a (k)) = V na{k) (p a (k))e-hWplW (21) 

with V na (k){pa(k)) indicating a polynomial of order n a (k). The only possible value for ^(k) 
that is positive and allows for a normalizable eigenfunction is given by 

7(^ = ^(^-/3^)) (22) 

with 

A(k) = (3 2 e 2 (k)+A . (23) 

In Figured] we plot the function ^(k) showing that as a function of f3e{k) the function is 
monotonically decreasing. This implies that the distribution Ti na {k){Pa{k)) is broader as /3 
increases and therefore the average energy (U({p(k)})) will increase with (3. In the limit 
/9e(/c) <C 1 the function 'y(k) given by Eq. (|2"2"|) converges to its maximum 

j(k) = ^ + <W*0) • (24) 
In the limit /3e(A:) ^> 1 instead we find the asymptotic behavior 

^ k)= ^[-^) + °(^m)\ (25) 
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FIG. 1: The function ^(k) denned in Eq. (I24p as a function of (3e(k), where e(k) has been defined 
in Eq. (|17p and j3 is the inverse temperature. 

The eigenvalue A^ of the eigenvalue problem Eq. (l2Tj) associated to the eigenfunction 
7Tfi({p(k)}) determined in Eq. (120]) is given by 

= II [1 - PmcMkt R{k)+niik)+1] (26) 

k 

Therefore the eigenvalues A^ only depend on the variables N k = nn(k) + rii(k). We can 
write these eigenvalues as 

A**) = e-^(W» (27) 

with 

E({N(k)}) = EM^k) 

k 

= ~\ £ + !] M 1 - Prnc^k)]. (28) 

Since N(k) = nn(k) + ni(k) each value E k (N k ) has a degeneracy gk(N k ) = (N k + 1). In 
Figure[2]we plot Ek(N k ) for iV^ = 0, 1, 2, 3, 4 versus the inverse temperature f3 while in Figure 
[3] we plot the same functions versus the energy e(k). From these figures it is clear that the 
levels Ek(Nk) of our model do not cross as a function of the inverse temperature (3. This is 
one of the major differences to traditional quasi-species models in which the level crossing is 
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associated with the error-threshold phase transition. While in biological evolution driven by 
a single-peak landscape, the system of biologically replicating entities is not able to localize 
on the fitness peak if the mutation rate is too high, in our Gaussian model of evolution 
of field units, this phase transition does not occur and the maximal eigenvalue associated 
to the fundamental state is an analytic function of /3. This is dependent on the particular 
quadratic shape of the energy function we have assumed, and other choices of the energy 
functional U({p(x)}) might lead to a phase transition. In the limit 0e(k) <C 1 the function 
E(Nf-) given by Eq. (1281) takes a relativistic form since using the definition of 'j(k) given by 
Eq. d22D and A(k) given by Eq. (F23]) we get 

E k (N k ) = ^/(hkc) 2 + (mc 2 y[N(k) + 1)] (29) 
with N(k) = nn(k) + nj(k). In the limit (3e(k) > 1 we obtain instead 

E (i N *}) = \ EM + !) ln [MQ/y/2\ . (30) 

Let us complete the solution of the stochastic evolution described by Eq. ( flOl) . First we 
decompose the function P({p(x)},r) in the basis ^n({p(k)}) given by Eq. ( 1201) and we get 

P({p(k)},r) = Y, c n( k )(r)n Ti ({p(k)}). (31) 

n(fe) 



Since the dynamical Eq. ( flOl can be easily linearized as shown in |2j, the solution of Eq. 
( TTUj) is immediate. Using the definition of the eigenvalue = exp[— (3E({N k })] given by 
Eq. (EID we get 



~n(fc)(r) = exp 



e -mm}) T _ z ] cm(0) _ (32) 



The function Z T appearing in Eq. (132]) and the function Z T defined in Eq. ( TTT]) are given 



by 



(33) 



Z T = J2e- mMk)}) c m (t) 

n(k) 



-2pE{{n(k)}) 



where (. . .) is the average over the ensemble of scalar field units. The evolutionary process 
described by Eq. (flO|) is an off-equilibrium and dissipative. In fact using Eq. ( |33|) and the 
dynamical solution for ca(k) given by Eq. ( 1321) we get 



dZ T 

~dT 



-20E({n(k)}) 



-(3E({n(k)}) 



(34) 
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FIG. 2: The first five levels Ek(Nk) as a function of the inverse temperature f3. The function is 
plotted for e{k) = 1 which is defined in Eq. ()17|) . 

Therefore, asymptotically in time the ensemble of scalar fields concentrates on the maximal 
eigenvalue of the evolutionary dynamics, similarly to what happens in biological evolution 
under the fundamental theorem of natural selection In fact the ratio between the time 
dependent coefficients c^ir) , Cft 2 {r) follows the relation 



chM 

CnM 



cxp 



n e -^({ni(fc)}) _ e -/3B({n 2 (fc)})j r j _ ^5) 



From Eq. f l32|) it is easy to show that the eigenfunction corresponding to the largest eigen- 
value Xft(k) an d lower value of E({Nk}) is dominating at large time scales. Interestingly, 
if the probability distribution P({p(k)},0) has a non-zero contribution c H= ^(0) > of the 
non-degenerate fundamental state {iV fc } = (0, . . . , 0), the probability distribution asymp- 
totically in time converges to the fundamental eigenfunction 



P({p(A;)},r)^n II MPa(k)) (36) 

k a=R,I 



corresponding to the value 



E({N k = 0}) = ~ £ ln[l - /3mc 2 7(A;)] . (37) 

' k 

If we impose a cutoff K to the momentum k (i.e. k < K ) we get that in the limit (3e(K) 1, 
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FIG. 3: The first five levels Ek{N k ) as a function of e(k) defined in Eq. (|17[) . The function is 
plotted for fixed inverse temperature 13 = 1. 

E({N k = 0}) takes the form 



E({N k = 0}) = J2 V ( hkc ) 2 + ( mc2 ) 2 ' ( 38 ) 

k 

In the opposite limit in which fimc 2 3> 1, E{{N k = 0}) takes the form 

E({N k = 0}) = §£ln (fie(k)/V2) . (39) 

The fundamental state is characterized by the eigenfunction in Eq. ( 136|) with iio(p a (k)) 
given by a Gaussian as described by Eq. (]2ip . Therefore, using the definition of 7(fe) given 
by Eq. (P2J we find 

2\ 2mc 2 1 

(Pj = 7W = y/A(k)-l3e(ky { } 

with a = R, I. Finally, asymptotically in time, the average energy (U({p(k)})) of the 
field units is given by the average of U({p(k)}) defined in Eq. ffT6]) on the fundamental 
eigenfunction given by Eq. (l36lh Therefore using Eq. ( 1401) we find 

(U({p(k)})) = V . 2e{k) . (41) 

In the limit /3e(K) < 1 we find that the average energy (U({p(k)})) is linear in e(k), 



(U({p(k)})) = J(Kkc)* + (mc^. (42) 
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Therefore in this limit, we find that, asymptotically in time, the average energy of the field 
units (U({p(k)})) is linear in e(k) acquiring in this way an apparently relativistic expression. 
In the limit f3mc 2 ^> 1 instead, we find that (U({p(k)})) is quadratic in e(k), 



M{pW}))=/3E[(^ c ) 2 + (^ 2 ) 2 ]. (43) 



k 

This shows that the stochastic quantization we have explored allows both for relativistic and 
non-relativistic expression of the average energy (U({p(k)})). 

V. CONCLUSIONS 

Erwin Schrodinger [3^] proposed that the complexity of living systems is related toquan- 



tum mechanics. Since then this idea has fascinated both physicists and biologists [3l|, |32 |. 
From the publication of What is life? however more than sixty years have passed and we 
still lack a solid scientific basis for Schrodinger's proposal. However the question Is life 



physics? 33] is attracting the interest of an increasing number of scientists. The connec- 
tion between the stochastic dynamics, which lies at the heart of biological evolution, and 
quantum-like behavior opens the way to apply the powerful techniques of the formalism of 
quantum mechanics within the theory of evolution. Recently the literature on this topic is 
acquiring a certain momentum and path-integrals, Schrodinger equation in imaginary time, 
Fock formalism and Bose-Einstein condensation are techniques which have been applied to 



study the quasi-species equation and stochastic effects in biological evolution [4Hl6|. These 
techniques do not solve all the complex aspects of biological evolution, i.e. biological evo- 
lution is not exclusively described by quantum formalism, but these techniques have been 
demonstrated to be extremely useful for solving specific biological questions. 
Many interesting works relate stochastic dynamics and quantum mechanics using stochastic 



quantization techniques 19M2 lL l25j. In this paper we have investigated at the same time the 
stochastic nature of biological evolution and the relation of this stochastic dynamics with 
stochastic quantization. The quasi-species equation describes the evolution of the proba- 
bility that a random individual in a population carries a given genome. Here we map the 
quasi-species equation for individuals of a self-reproducing population to an ensemble of 
scalar field elementary units undergoing a creation and annihilation process. In this map- 
ping the field units correspond to individuals and the field p(x) to their genome. The selective 
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pressure instead maps to an "inverse temperature" /3 that parametrize the evolution of the 
field units. We have solved the specific case in which the role of the Fisher fitness is played 
by a quadratic energy potential of the field and the mutations by Gaussian noise on the 
field values. The ensemble of field units described by the quasi-species equation relaxes to 
the fundamental state, describing an intrinsically dissipative dynamics. We have presented 
interesting relations between this stochastic evolution and quantum-like behavior with the 
emergence of discrete eigenvalues for the particle probability distributions. For a quadratic 
dispersion relation the mean energy (U) of the system changes as a function of the inverse 
temperature /3. For small values of /3 the average energy (U) takes a relativistic form, for 
large values of /3, the average energy (U) takes a classical form. 

In future works we plan to study how this mean-field picture described by the quasi-species 
equation is modified by stochastic effects. 

G.B. acknowledges interesting discussions with H. Goldberg, G. Jona Lasinio, M. Kardar, 
P. Nath, M. A. Nowak, and S. Redner. 
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